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Results of a study to validate the Fast Scattering Code (FSC) as a duct noise predictor, includ- 
ing the effects of curvature, finite impedance on the walls, and uniform background flow, are 
presented in this paper. Infinite duct theory was used to generate the modal content of the 
sound propagating within the duct. Liner effects were incorporated via a sound absorbing 
boundary condition on the scattering surfaces. Simulations for a rectangular duct of constant 
cross-sectional area have been compared to analytical solutions and experimental data. Com- 
parisons with analytical results indicate that the code can properly calculate a given domi- 
nant mode for hardwall surfaces. Simulated acoustic behavior in the presence of lined walls 
(using hardwall duct modes as incident sound) is consistent with expected trends. Duct curva- 
ture was found to enhance weaker modes and reduce pressure amplitude. Agreement 
between simulated and experimental results for a straight duct with hard walls (no flow) was 
excellent. 
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Nomenclature 

Wave amplitude; also, complex admittance, A = 1/Z 

Wave amplitude 

Waveguide dimension 

Mach number 

Specific acoustic resistance 

Region inside rectangular waveguide 

Specific acoustic reactance 

Specific acoustic impedance, Z = R+Xi 

Speed of sound inside duct 

Excitation frequency, Hz 

Imaginary unit, 

Monopole index 
Wave number 
Pressure 
Time 

Cartesian coordinates 

Kronecker delta 

Air density inside duct 

sgn (z-z J Q ) , a notational convenience 

Excitation frequency, 2nf 

General coordinates 

Location of monopole source 
Order of mode in x direction 
Order of mode in y direction 
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x,y,z = Cartesian coordinate directions 

a = Monopole strength 

°o = Freestream quantity 

Superscripts 

' = Acoustic variable 

—> = Vector 

I - Introduction 

An accurate determination of sound as it propagates within waveguides of arbitrary geometry is of paramount 
importance to the development of noise reduction concepts for turbofan bypass ducts. Of interest are the effects asso- 
ciated with changes in duct shape, wall acoustic treatment, and the presence of components such as splitters and out- 
let guide vanes. 

An experimental study of noise abatement concepts for bypass ducts is costly and time consuming. Taking 
advantage of a recently developed computational method for predicting the scattered acoustic field generated by the 
interaction of a known source and an arbitrary collection of surfaces 1 , a multitude of configurations can be explored 
in a timely manner, giving the engineer the ability to select only the best candidates for experimental evaluation. The 
prediction methodology, which uses the equivalent source method (ESM) to solve an exterior Helmholtz equation 
boundary value problem (BVP), has been implemented into a computer program called the Fast Scattering Code 
(FSC). The code uses pre-defmed component geometries and operating conditions to simulate the time-harmonic 
scattered acoustic pressure and velocity fields at user specified locations. Details on the methodology and program 
usage are given in references 1 through 3. 

The main objective of this investigation is to validate the FSC as a duct noise predictor, including the effects of 
duct curvature, acoustic treatment, and background flow. To do so, simulations for a duct with a constant, rectangular 
cross-section will be compared to analytical results and experimental data obtained with the NASA Langley Research 
Center (LaRC) Curved Duct Test Rig (CDTR) 4 . For the present study, the coordinate system follows FSC conven- 
tion, that is, x increases “up”, y increases spanwise, and z increases upstream. The dimensions of the duct are similar 
to those of the CDTR (X x , Ly, L z = 0.381 m, 0.1524 m, 3.0 m). The length of the test section (where the liner speci- 
mens are placed) is 0.8128 m. For simplicity, the origin lies on the duct axis, at the center of the test section. Har- 
monic waves of the form e mt are assumed. 

II - Validation with Analytical Results 


A. Infinite Duct Theory 

In general, the simulation of a scattered acoustic field using the FSC requires discretization of the surface(s) 
according to excitation frequency and surface area, and the generation of a suitable sound source. For constant area 
duct applications, the former requirement is easily performed using grid generation software; the latter can be simu- 
lated using combinations of simple monopole sources to produce duct modes of arbitrary order, amplitude, and phase. 
Consider a waveguide of infinite length, centered at the origin: 
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By judicious choice of (x 0 , yg), certain modes can be silenced by making 4 mn (x 0 ,y 0 ) = 0 (m and n not both 
zero). For example, to silence the n-th mode,y 0 can be chosen so that 
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For J monopoles with strengths a,-, the acoustic pressure can be expressed as: 
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To silence multiple modes in regions of space for which z < min z{ 
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For example, two sources can be used to silence the (0,0) mode: 
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Various dominant duct modes, excitation frequencies, background flow speeds, and wall impedances were used 
in the present investigation. The modes, frequencies, and flow settings were chosen from a subset of experimental 
data obtained using the CDTR with a straight test section; the impedances were educed from measurements taken at 
the NASA LaRC Grazing Incidence Tube (GIT) facility 6 for a single degree of freedom, locally reacting liner built 
in-house. The liner, denoted as L01 7 , has a honeycomb core 1.5 inches deep covered with an 8.7% open perforate. 
The parameter values are summarized in Table 1 . 


Table 1 . - Parameters used in validation study. 


Excitation frequency 

1500 Hz 

2000 Hz 

Duct mode 

(0, 0) 

(2,0) 

(2, 1) 

Normalized Liner 
impedance, = 0.0 

°° + O.Oi 

0.23 -0.14i 

°° + O.Oi 

0.23 -0.14i 

°° + O.Oi 

0.12 + 0.35i 

Normalized Liner 
impedance, = 0.275 

°o + O.Oi 

0.76 - 0.43i 

oo + O.Oi 

0.76 - 0.43i 

go + O.Oi 

0.73 + 0.13i 


At an excitation frequency of 1500 Hz, several modes are cut on: (0, 0), (1, 0), (2, 0), (0, 1), (1, 1), (3, 0), and (2, 
1). The n= 1 modes, (0, 1), (1, 1) and (2, 1), are silenced by placing one source at = 0-0; the odd m modes, (1, 0) 
and (3, 0), are silenced by placing the source at x 0 = 0.0. The remaining modes, (0, 0) and (2, 0), are alternately 
silenced using equation (10). At an excitation frequency of 2000 Hz, two additional modes are cut on, (3, 1) and (4, 
0). Four monopoles can be used to isolate the (2, 1) mode. Placing the sources at y 0 * 0 ensures generation of the n = 
1 modes; positioning one pair of sources at x 0 = ±3L x /8 silences the (1, 0), (3, 0), (1, 1), and (4, 0) modes. The addi- 
tion of a second pair of sources at the same (x 0 , yg) but different z 0 silences the (2, 0) mode. Modal decomposition 
indicates that, of the remaining modes, (0, 0), (0, 1), and (2, 1), the amplitude of the (2, 1) mode is approximately 5 
times larger than that of the next strongest one. 
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B. Straight Duct 

Mode propagation patterns at the walls and bisecting plane of a hard-wall, straight duct for excitation frequencies 
of 1500 Hz and 2000 Hz, in the absence of flow, are shown in Figure land 2 for dominant (0, 0), (2, 0), and (2, 1) 
modes. The axial placement of the acoustic sources, z J and z^ , corresponds to the locations of the two driver planes 
in the CDTR rig. The test section inlet (in) and outlet (out), indicated in the figures, are at z = 0.4064 m and -0.4064 
m, respectively. Although reflections from the duct openings may be present, note from Figure la that the FSC is able 
to simulate the proper behavior of a dominant, aft travelling plane wave. Higher order modes are correctly simulated 
as well. Figure lb indicates that a clearly dominant (2, 0) mode is established immediately aft of the monopole source 
combination. The presence of several modes can be observed in Figure 2. As indicated by infinite duct theory, the (2, 
1) mode appears dominant on the walls of the duct, and a plane wave of amplitude approximately 1/5 that of the (2, 1) 
mode dominates the propagation behavior at the bisecting plane, which coincides with a node of the (2, 1) mode. 
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Figure 1. - Mode propagation within a hardwall, straight duct; /= 1500 Hz, M ^ = 0.0. 
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Figure 2. - Mode propagation within a hardwall, straight duct ;/= 2000 Hz, M ^ = 0.0, dominant (2, 1) mode. 
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The experimental rig is equipped with a flow system capable of delivering air at speeds typical of aft fan bypass 
duct flows. The effects of background flow (M^ = 0.275) on noise propagation inside the straight duct are depicted in 
Figures 3 and 4. In the laboratory-fixed frame of reference, air flows from right to left (i.e., in the -z coordinate direc- 
tion). For the same source locations as above, note from Figures 3a, 3b, and 4 that the main effects of flow appear to 
be 1) a reduction in the pressure amplitudes inside the duct, 2) a decrease in wave length upstream of the sources, and 
3) an increase in wavelength downstream of the sources. A point monopole source at rest emits an omnidirectional 
sound field. If the source is in motion, the radiated sound field is distorted so that pressure perturbations forward of 
the source are more intense 8 (with reduced wavelengths) than those aft of it (with increased wavelengths), as 
observed in the predictions. 
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(b) Dominant (2, 0) mode 


Figure 3. - Mode propagation within a straight duct with hard walls; /= 1500 Hz, 
M ^ = 0.275 (air flows in the -z direction). 
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Figure 4. - Mode propagation within a straight duct with hard walls; /= 2000 Hz, 
M ^ = 0.275 (air flows in the -z direction), dominant (2, 1) mode. 
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The sound absorbing boundary condition imposed on the scattering surfaces 3 requires the input of dimensional 
admittance, A. To simulate hard walls, A = 0.0 + 0.0/. To simulate acoustically treated walls, the normalized imped- 
ances of Table 1 are converted to admittances using environmental variables (air density, p, and sound speed, c, inside 
the duct) measured during CDTR tests. For the present exercise, the incident acoustic field is in the form of hardwall 
duct modes. 

The effects of liner treatment for the entire side walls of the test section (the top and bottom walls remain hard) 
are shown in Figures 5 through 7 for M ^ = 0.0. For a given frequency, the peak attenuation of a liner is proportional 
to the length and surface area of the treatment 9,10 . It has also been determined 9 that the initial portion of the treat- 
ment, in the vicinity of the test section entrance, is generally more effective because the higher order modes are atten- 
uated more rapidly than the remaining - predominantly plane-wave - propagating modes. Note from Figure 5 that, for 
a dominant (0, 0) mode, the effect of the liner increases gradually from the test section inlet plane (z = 0.4064 m), 
peaks at approximately z = 0.0 m, and remains more or less constant after that. The sound attenuation through the 
treated test section appears to be on the order of 30 dB. Similar behavior is observed for a dominant (2, 0) mode (Fig- 
ure 6), with a slightly larger peak attenuation (35 dB). When the (2, 1) mode is dominant in the duct, the liner seems 
to abate all modes except for the plane wave, as clearly seen in the acoustic pressure contours of Figure 7. This results 
in a relatively modest sound reduction through the test section, approximately 5 dB. 



Figure 5. - Mode propagation within a straight duct with partially lined side walls, Z= (0.228 - 0.137/) (pc); 

/= 1500 Hz, = 0.0, dominant (0, 0) mode. 
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Figure 6. - Mode propagation within a straight duct with partially lined side walls, Z = (0.228 - 0.137/) (pc); 

/= 1500 Hz, = 0.0, dominant (2, 0) mode. 




Figure 7. - Mode propagation within a straight duct with partially lined side walls, Z= (0.1 15 + 0.347/) (pc); 

/= 2000 Hz, = 0.0, dominant (2, 1) mode. 
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The effects of a uniform background flow on liner performance are presented in Figures 8 through 1 1 . For a dom- 
inant plane wave at /= 1500 Hz (Figure 8), the introduction of flow appears to enhance the effectiveness of the liner 
in the region near the test section entrance. Most likely, this is a result of the changes in wavelength and/or pressure 
amplitude that occur in the presence of flow. The peak attenuation through the test section did not seem to be affected, 
remaining at about 30 dB. No appreciable effects on liner behavior due to flow for a dominant (2, 0) mode at /= 1500 
Hz can be observed (Figure 9). The absence of noticeable changes may be caused by the imposition of a non-opti- 
mum liner impedance on the test section walls - the same value educed for the dominant plane wave was used for the 
(2, 0) mode as well, since the experimental data set did not contain information for this mode/frequency combination. 

Flow effects on liner performance for a dominant (2, 1) mode at /= 2000 Hz are depicted in Figures 10 and 11. 
When the surface impedance does not include the effects of flow (Figure 10), the change in acoustic behavior inside 
the duct is similar to that for a dominant plane wave: the liner is slightly more effective near the test section inlet. 
When the effects of flow are included in the impedance (Figure 11), liner performance is noticeably improved, result- 
ing in further attenuation of the downstream propagating plane wave. This effect is apparent in the acoustic pressure 
contours of Figure 11. Note also from the SPL contours that the attenuation through the test section increases to 
approximately 20 dB. It has been determined 7 that the presence of flow alters surface impedance by increasing resis- 
tance, R , and shifting reactance, X , to higher frequencies. The magnitude of the shift depends on the location of a par- 
ticular excitation frequency on the X(f) curve for the L01 liner: 1500 Hz may fall on a region of slowly varying 
reactance, that is, near the resonant frequency of the liner, whereas 2000 Hz may be within a region of rapidly chang- 
ing reactance. The addition of flow also shifts peak attenuation to higher frequencies, which may result in a consider- 
able reduction in plane wave amplitude, as seen in Figure 11. 
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Figure 8. - Mode propagation within a straight duct with partially lined side walls, Z = (0.761 - 0.427/) (pc); 
/= 1500 Hz, M ^ = 0.275 (air flows in the -z direction), dominant (0, 0) mode. 
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Figure 9. - Mode propagation within a straight duct with partially lined side walls, Z= (0.761 - 0.427/) (pc); 
/= 1500 Hz, M ^ = 0.275 (air flows in the -z direction), dominant (2, 0) mode. 
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Figure 10. - Mode propagation within a straight duct with partially lined side walls, Z= (0.115 + 0.347/) (pc); 
/= 2000 Hz, M ^ = 0.275 (air flows in the -z direction), dominant (2, 1) mode. 
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Figure 11. - Mode propagation within a straight duct with partially lined side walls, Z= (0.730 + 0.130/) (pc); 
/= 2000 Hz, M ^ = 0.275 (air flows in the -z direction), dominant (2, 1) mode. 


C. Duct with Curved Test Section 


A primary goal of the CDTR project is the evaluation of duct curvature effects on noise propagation and liner 
performance. This is accomplished by modifying the test section so that the outlet plane is offset a distance of 0.0762 
m or 0.1524 m , that is, 50% or 100% of the duct width, L y . Simulations for a hardwall curved test section with an off- 
set of 1.0 L y are presented in Figure 12 for a dominant (2, 0) mode at an excitation frequency of 1500 Hz, = 0.0. 
Note from the figure that interactions with the curved surfaces cause a slight decrease in wave amplitude aft of the 
test section. No appreciable change in liner behavior as a result of duct curvature can be observed (Figure 13). 



Figure 12. - Mode propagation within the CDTR test section with hard walls, duct offset =1.0 
L y ; /= 1500 Hz, M ^ = 0.0, dominant (2, 0) mode. 
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Figure 13. - Mode propagation within the CDTR test section with lined walls, Z = (0.228 - 0.137z) (pc), 
duct offset =1.0 L y \ f= 1500 Hz, M ^ = 0.0, dominant (2, 0) mode. 


The effects of a more aggressively curved test section on mode propagation are depicted in Figures 14 and 15 for 
/= 1500 Hz, M ^ = 0.0. The curvature follows a full sine wave of amplitude 0.5 L y . Note from Fig. 14 that the charac- 
teristics of a dominant plane wave propagating inside a hardwall duct are considerably altered inside the test section: 
secondary modes become apparent, causing an increase in wavelength at the walls and a decrease at the center. A 
dominant plane wave of reduced amplitude is re-established aft of the test section. Curvature causes a slight enhance- 
ment of the weaker modes inside the test section when a dominant (2, 0) mode is generated, as seen in Figure 15. 
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Figure 14. - Mode propagation within a sinusoidal test section with hard walls;/= 1500 Hz, 

= 0.0, dominant (0, 0) mode. 
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Figure 15. - Mode propagation within a sinusoidal test section with hard walls;/= 1500 Hz, 

= 0.0, dominant (2, 0) mode. 
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The characteristics of a wave traveling down the axis of a hardwall duct, as influenced by test section curvature, 
are depicted in Figure 16 for a dominant (2, 0) mod e,/= 1500 Hz, M ^ = 0.0. At this location within the duct, the 
salient features are a reduction in wave amplitude with increasing curvature, accompanied by changes in wavelength. 
Thus, it appears that duct curvature attenuates modal noise propagating within a waveguide of rectangular cross-sec- 
tion. Liner behavior does not seem to be affected by increased curvature of the test section (see Figures 17 and 18). 



Figure 16. - Effect of test section curvature on a dominant (2, 0) mode, hardwall duct; /= 1500 Hz, M ^ = 0.0. 



Re (p). Pa 


Re (p). Pa 


Figure 17. - Mode propagation within a sinusoidal test section with lined side walls, 
Z= (0.228-0.1370 (pc);/= 1500 Hz, = 0.0, dominant (0, 0) mode. 
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Figure 18. - Mode propagation within a sinusoidal test section with lined side walls, 
Z= (0.228-0.137/) (pc);/= 1500 Hz, = 0.0, dominant (2, 0) mode. 
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III. Validation with Experimental Results 


The exterior Helmholtz BVP, which is solved using a collection of point multipoles located inside the scattering 
body (equivalent sources), can be transformed into an interior Helmholtz BVP by placing the equivalent sources out- 
side the scattering surface. Doing so can substantially reduce the size of the BVP, since only the inner surfaces of the 
duct need to be considered. A Dirichlet boundary condition for initialization of acoustic pressure at the duct openings 
has been implemented into the FSC. 

For this investigation, the incident sound field was input in the form of total acoustic pressure at the inlet (z = 
0.4064 m) and outlet (z = -0.4064 m) planes of the test section. The pressures at these locations were synthesized 
from experimental measurements obtained with the CDTR for a target amplitude of 60 Pa (130 dB) and a mode 
amplitude ratio (MAR)^ > 10 dB. The size of the input planes, 12 x 30, is sufficient for proper data resolution up to 
2400 Hz (current frequency limit for the rig). 

The results presented here correspond to a test section with hard walls, without flow. Although experimental data 
were available for lined side walls, the impedances educed with the GIT are incompatible with the formulation of the 
interior BVP used for this study. A proper comparison between simulated and experimental results for ducts with soft 
walls can be performed only if the impedances are educed with the FSC using the given synthesized pressures at the 
duct openings. Such an endeavor is beyond the scope of this paper. 

Simulated acoustic pressures at the walls and bisecting plane of an unlined straight test section for a dominant (0, 
0) mode at /= 1500 Hz, in the absence of flow, are depicted in Figure 19. Note that the code is able to calculate the 
proper modal behavior inside the duct using the Dirichlet boundary condition. Comparisons between calculated and 
measured acoustic pressures at the duct openings are presented in Figure 20. The agreement is excellent. 

Simulations for a dominant (2, 0) mode at 1500 Hz, without flow, are presented in Figures 21 and 22. Note from 
the pressures at the surfaces and bisecting plane (Figure 21) that the correct mode propagates through the test section. 
Agreement between calculated and experimental pressures at the duct openings (Figure 22) is excellent. 



Figure 19. - Wave propagation within straight CDTR test section with hard walls ;/= 1500 Hz, 

= 0.0, dominant (0, 0) mode. 


{. MAR is the ratio of the target mode’s amplitude to the amplitude of the next strongest mode, expressed in 
dB. The target mode is considered dominant if MAR > 10 dB. 
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Figure 20. - Comparison of total acoustic pressure at the test section openings, hard walls; 
/= 1500 Hz, = 0.0, dominant (0, 0) mode. 
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Figure 21 . - Wave propagation within straight CDTR test section with hard walls; /= 1500 Hz, 

= 0.0, dominant (2, 0) mode. 
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Figure 22. - Comparison of total acoustic pressure at the test section openings, hard walls; 
/= 1500 Hz, = 0.0, dominant (2, 0) mode. 
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IV. Concluding Remarks 


A study aimed at validating the Fast Scattering Code as a duct noise predictor has been undertaken. The effects 
of duct curvature, acoustic treatment, and uniform background flow have been included. Hardwall duct modes were 
used as the incident acoustic field for exterior BVP calculations. Comparisons between analytical and simulated 
results for a selected set of modes and frequencies indicate that the FSC is capable of generating the correct modal 
behavior inside the duct. Liner effects were incorporated through the use of a sound absorbing boundary condition at 
the scattering surfaces. The main effects of liner technology - changes in modal content and a reduction in acoustic 
pressure amplitude - are apparent in the simulations. The addition of curvature to the test section resulted in a local- 
ized amplification of weaker modes, and a reduction in amplitude for aft propagating waves. A Dirichlet boundary 
condition for use with interior BVPs has been added to the code. Comparisons between simulated and experimental 
results for a straight, hardwall test section without flow are excellent. 
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